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Abstract 

Some new aspects of axially symmetric spacetimes are discussed. 
These results open the door for future interplay between analytical 
and numerical studies. The new developments are based on the role 
of the total mass in axial symmetry. Finally, a list of relevant open 
problems is presented. These problems can be hopefully solved with 
an interaction between numerical and analytical insights. 

1 Introduction 

In any physical theory, the presence of a symmetry reduces the degrees of 
freedom of the equations and hence it simplifies considerable the analysis. To 
study an isolated systems, the simplest models are the spherically symmetric 
ones. However, it is well known that for vacuum Einstein equations, due to 
Birkhoff's theorem, spherical symmetric spacetimes has no dynamics. The 
next possible model with symmetries in vacuum are axially symmetric space- 
times. It has been proved in [1] that no additional symmetry can be imposed 
to the spacetime if we want to keep the gravitational radiation and a com- 
plete null infinity. This result single out axially symmetric spacetimes as the 
only models for vacuum isolated, dynamical, system with symmetries. From 
this point of view, axially symmetric gravitational waves are the simplest 
possible waves emitted by isolated sources. 
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There exists many relevant physical models one can study in axial sym- 
metry. In particular, for vacuum and in the strong field regime, we list the 
following 

• Head-on collisions of two black holes. 

• Rotating, non-stationary, black holes. 

• Formation of black holes: weak cosmic censorship. 

• Critical collapse of gravitational waves. 

What are the difficulties of axial symmetry? To take advantage of the 
symmetry an adapted coordinate should be used. However, the reduced 
equations are formally singular at the axis. This singular behavior introduce 
a major difficulty in the analysis. In fact, it can be argued that this singular 
behavior near the axis is so complicated that the axially symmetric case is 
as hard as the full general case. 

What are the advantages of axial symmetry? The first obvious advan- 
tage for numerical computations is that axially symmetric spacetimes are less 
computationally expensive: only two effective spatial dimensions. Also, the 
number of equations and variables is reduced. The second advantage, which 
is also well known, is the conservation of angular momentum. The angular 
momentum is a quasilocal conserved quantity in axial symmetry (Komar in- 
tegral of the Killing vector). That is, axially symmetric gravitational waves 
do not carry angular momentum. In particular no Penrose process and no 
superradiant scattering can occur for axially symmetric vacuum spacetimes. 
This represents an important simplification in the dynamics. The third ad- 
vantage is the mass integral formula. The total mass can be written as a 
positive definite volume integral, as we will see. This is the main new ingre- 
dient that we would like to discuss in the following. 

The purpose of this review is to summarize some new results for axi- 
ally symmetric spacetimes with emphasis in the interplay between numerical 
and analytical studies. The plan of the article is the following. In section 
[2] we review the axially symmetric equations. The purpose of this section 
is to present the relevant maximal-isothermal gauge and the corresponding 
mass formula, which is valid only in this gauge. This formula represents 
the main new motivation to study axial symmetry. In section [3] we dis- 
cuss recent results, both numerical and analytical which are based on the 
maximal- isothermal gauge and the mass formula. Finally, in section H] we 
present relevant open problems. 
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2 Axial Symmetry 



Consider a vacuum solution of Einstein's equations, i.e., a four dimensional 
manifold M with metric g^ (with signature ( — h ++)) such that the corre- 
sponding Ricci tensor vanishes 

{4) n, u = o. (i) 

Suppose, in addition, that the metric g^ u admits a Killing field rf ', that is rf 
satisfies the equation 

V(^) = 0, (2) 

where is the connection with respect to g^. Greek indices //, u, • • ■ denote 
four dimensional indices. 

We define the square of the norm and the twist of respectively, by 

V = rfrfg^, w„ = e^^V V- (3) 
Using the field Eq. (pQ) it is possible to prove that 

V^j = 0, (4) 
and hence is locally the gradient of a scalar field u 

Up = Vf.co. (5) 

Let M denote the collection of all trajectories of rf , and assume that it 
is a differential 3-manifold. We define the metric h^ v on M by 

W\iv = h^ u + r]^rj u . (6) 

The vacuum field equations ([1]) can be written in the following form on M 

n v =-(V a V V aV -V a uV a u), (7) 
V 
2 

Uuj = -V a uVaV, (8) 
V 

{3) Rab = T^iVaV^bV + V a uV b 0j), (9) 

where V a and ^R a b are the connexion and the Ricci tensor of h a b, we have 
defined □ = V a V a and Latin indices a,b . . . denote three dimensional indices 
on J\f. 
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Up to this point, the only assumption we have made is that the spacetime 
admits a Killing vector field rf and that 77^ is not null, otherwise the metric 
h a b is not defined. If the Killing field is timelike (77 < 0) then the metric h a b 
is Riemannian and the equations dZD - © ar e the stationary Einstein vacuum 
equations. On the other hand, when the Killing vector is spacelike (77 > 0), 
the metric h ab is a is a 3-dimensional Lorenzian metric (we chose the signature 
( — h +)). In axial symmetry, the Killing vector rf is spacelike and its norm 
vanishes at the axis of symmetry. Hence, the equations are formally singular 
at the axis. This singular behavior at the axis represents the main difficulty 
to handle these equations. 

In the Lorenzian case, Eq. (Q has the form of Einstein equations in 
three dimensions coupled with effective matter sources produced by 77 and u. 
The effective matter Eqs. (JTJ)— ([HD imply that the energy-momentum tensor 
defined in terms of 77 and oj by 

Tab = 2-^(V a 77V 6 7? + VauVbUj) - -^h ab {V c rfJ c V + V c uV c u), (10) 

is divergence free, i.e. V a T ab = 0. 

Eqs. dZD-flSJ) are purely geometric equations with respect to the metric 
h a b. The essential point is that there exist no dynamical degrees of freedom in 
3-dimensional gravity (the Weyl tensor vanishes) and hence all the dynamics 
is produced by the effective matter sources determined by (t),lo). In other 
words, the dynamics of vacuum axially symmetric gravitational waves has 
the behavior of matter moving in a lower dimensional space. Hence, it can 
be expected that the total mass of the waves can be computed as an "integral 
on the matter sources" . As we will see in the following, this is precisely what 
happens and it leads to the mass formula in axial symmetry. 

2.1 2+1 decomposition 

In order to formulate an initial value problem, we will perform an standard 
2 + 1 decomposition of Eqs. ((ZD - ©- Note that this is completely analogous 
to the 3 + 1 decomposition of Einstein equations, in fact all the formulas are 
formally identical because the dimension does not appear explicitly in them 
(see, for example, [28], [27]). 

Consider a foliation of spacelike, 2-dimensional slices S of the metric h ab . 
Let t be an associated time function and let n a be the unit normal vector 
orthogonal to S with respect to the metric /i &. The intrinsic metric on S is 
denoted by q a b and is given by 

hab = -n a n b + q ab . (11) 
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Define the density fx by 

H= 2 {3) R ab n a n b + ^R, (12) 

and the current J b by 

J b = -qln a ®R ca , (13) 

where ( 3 >R = < ^ > R a bh ab denotes the trace of ^'R a b- Then, using Eq. ([9]) we 
obtain 

V = ^ 2 W 2 + + \DV\ 2 + \Du\ 2 ) , (14) 

J a = {v'DaV + ui'D A u) , (15) 

where Da is the connexion with respect to qAB- The prime denotes directional 
derivative with respect to n a , that is 

r{ = n a V aV = - (d tV - (3 A D A ri) (16) 
a 

where a is the lapse and j3 A is the shift vector of the foliation. The indices 
A, B, ■ ■ ■ denote two dimensional indices on S. The constraints equations 
corresponding to are given by 

{2) R-x ab xab + x 2 = ^ (17) 

D A X ab - D BX = Jb, (18) 

where ^R is the Ricci scalar of qAB, Xab is the second fundamental form of 
S and x hs trace 

X = q AB XAB- (19) 
We use the following sign convention for the definition of xab 

Xab = -q c a S7 c nb = ~^£ n qab, (20) 

where £ denotes Lie derivative. The evolution equations are given by 

d t qAB = -2axAB + £/3qAB, (21) 

d t Xab = £?Xab ~ D A D B a + ar AB , (22) 

where 

TAB = XXAB + {2) RAB ~ {3) RAB ~ 2 X ACX%- (23) 
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and 

{3) Rab = ^(dAVdsV + d A ud B u). (24) 

The evolution equations fl2T]) - fl22|) and the constraint equations f|T7|) -f lT8|) 
constitute a complete 2 + 1 decomposition of the 3-dimensional Einstein Eq. 
([§]). It remains to decompose the effective matter Eqs. (JTJ)— (|HD - The result 
is the following 

D A a 1 

-£" + A 9 S + D a Tj—— + E' X = ±- W 2 ~ \Duj\ 2 ) , (25) 
D A a 2 

-u" + A q co + D A u + u' X =— (D A ujD A r] - co'r]') , (26) 

a r\ l v ' 

where A q is the Laplacian with respect to qAB, i-e. A q = D a Da and we have 
defined £ = log?]. 

The evolution equations (12"TT) - (|2"21 and fl23]) - (|2^]) together with the con- 
straint equations f JT7|) - f|T8|) represent the complete 2 + 1 decomposition of the 
geometrical equations (JTJ) — ([9]) in an arbitrary gauge. This set of equations 
presents an important feature which is peculiar of 2 + 1 dimensions. Namely, 
the evolution equations fl2"TT)- fl2"2l and the constraint equations f|T71) - ffT51) are, 
roughly speaking, equivalent in the following sense. With an appropriate 
gauge choice (in the next section we will present an important example) the 
constraint equations (TT71)- (|T8l) form an elliptic system that determines qAB 
an xab in terms of the effective matter sources r\ and u. Hence, the evolution 
equations (125]) -( 126]) together with the constraint equations ( TT7|) - (|T8|) consti- 
tute a complete system of equations. In particular, a solution of this system 
will automatically satisfy the other evolution equations (12 ip - (122 p . Alterna- 
tive, it is possible to solve the evolution equations ( 125]) - (T26"]) and f|2Tj) f|22|) . If 
the initial data satisfy the constraint equations at a given time, then the so- 
lution will satisfy the constraint equations for all times. Note that even when 
the gauge is fixed, we have different possibilities for constructing evolutions 
schemes. The first choice is called in the numerical literature "constrained 
evolution" and the second choice "free evolution" . It is also possible to mix 
them to obtain "mixed evolution" schemes (see the discussion in [29J). 



2.2 Gauge and mass formula 

In this section we describe the maximal-isothermal gauge. In particular we 
review the mass formula valid in this gauge (see [T3] for details). For the 
lapse, we impose the maximal condition on the 2-surfaces 

X = 0. (27) 
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Note that we are not imposing that the surfaces are maximal in the 3- 
dimensional picture as in [TJ]. The later condition is the one generally used 
(see, for example, j6] [30]), but the difference is only minor. In particular the 
mass formula is positive definite for both conditions. The one used here ap- 
pears to be natural with respect to the rescaled metric h a b- Eq. ( 127)) implies 
the following well known equation for the lapse 

A q a = a( X AB XAB + Vi), (28) 

where 

Ai 1 = < s >i^nV=^(V 2 + (29) 

The maximal gauge (127)) can be, of course, imposed in any dimensions and it 
is not related at all with axial symmetry. In contrast, the condition for the 
shift is peculiar for two space dimensions. The shift vector is fixed by the 
requirement that the intrinsic metric qAB has the following form 

Qab = e 2u S AB , (30) 

where Sab is a fixed (i.e. B^ab = 0) flat metric in two dimensions. Then, 
using (127)) . we obtain that the trace free part of (12Tj) is given by 

2axAB = {C q /3) A B, (31) 

where C q is the conformal Killing operator in two dimensions with respect to 
the metric qAB- Equation (13"T1) is an elliptic first order system of equations 
for /3 A . 

The elliptic Eqs. (1251) and (l3Tj) determine lapse and shift for the metric h a b 
and hence they fix completely the gauge freedom in Eqs. (JTJ) (J9j) . This gauge 
has associated a natural cylindrical coordinate system (t, p, z) for which the 
metric Sab is given 

S = dp 2 + dz 2 , (32) 

and the axis of symmetry is given by p = 0. The slices S are the half planes 
M 2 + . 

It is useful to introduce the following auxiliary functions. Instead of rj we 
will use the function a defined by 

V = pV. (33) 

Note that for Minkowski r\ = p 2 , that is, a measure the non-flat part of the 
norm rj. Also, instead of u it is convenient to work with the function q defined 
by 

u = log p + cr + q. (34) 
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The important property of q is that it vanished at the axis. This is a con- 
sequence of the regularity conditions of the metric at the axis (see [20] for 
details). 

The equations presented in section 12.11 can be explicitly written as partial 
differential equations in the coordinates (t,p,z) (see |20j). In particular the 
Hamiltonian constraint is given by 

^Aa + Aq = - 6 -, (35) 

where 

e = ^ W 2 + uj' 2 ) + |da| 2 + M + 2e~^x AB XAB. (36) 

The hat on \ AB means that the indices are moved with the flat metric Sab, 
d denotes partial derivatives with respect to p and z, A is the 2-dimensional 
flat Laplacian, namely 

A = d 2 p + d 2 z , (37) 

and ( 3 ^A is the 3-dimensional flat Laplacian acting on axially symmetric 
functions, that is 

(3) A = A + ^. ( 38 ) 
P 

The positive scalar e plays the role of the energy density of the gravita- 
tional waves in this gauge. The integral of e over the slice S is given by 

1 f 

m— — / epdpdz. (39) 

It is a non-trivial fact that this number m is precisely the total ADM mass 
of the spacetime. Moreover, the fall off properties of the solutions a and f3 
of gauge conditions ensure that the mass is conserved along the evolution, 
namely 

lm-0, (40) 

see [15] for details. 

The mass formula fl39|) together with the conservation law fj40|) represent 
a relevant property of this gauge which is not present in other ones. 

In any physical theory conserved quantities (in particular, conserved en- 
ergies) are very important to control the evolution of the system. However, 
in General Relativity, the conserved mass appears as a boundary integral and 
not as a volume integral (as, for example, in the wave equation). Hence it 
is not possible to relate the mass with any norm of the fields to control the 
evolution of them (for the wave equation the energy is precisely the norm 
of the wave). The mass formula fl39l) for axially symmetric systems in the 
maximal-isothermal gauge represents a remarkable exception. 
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3 Results 



3.1 Numerical results 

Axially symmetric spacetimes has been studied numerically since the very be- 
ginnings of numerical relativity, see chapter 10.4 in [2] and references therein 
for an historical perspective. A very good reference for early results is the 
review article [3]. In particular, in this article the maximal and isothermal 
gauges are discussed. Another important review article is [26J, where the 
symmetry reductions and 2 + 1 decomposition mentioned in section [2] was 
used for first time in numerical relativity. 

As we mentioned in the introduction, the difficulty introduced by the sin- 
gular behavior at the axis is severe and that made axially symmetric evolution 
particularly difficult to handle. As a consequence of this "axisymmetric codes 
were practically abandoned as soon as computers became powerful enough 
to be able to handle full three-dimensional simulations in the early 1990s." 
(the quote is from [2]). 

New insights in the study of the axial regularity problem in numerical 
relativity were presented in [23] • They introduced two main ingredients. The 
first one is that the behavior at the axis is treated as a singular boundary 
condition which are handled by introducing "ghost zones" outside the axis. 
The second one is that all the variables used are such that either they vanish 
at the axis or have vanishing derivative there, but not both. This cure the 
axis instabilities from the numerical point of view. 

Later in [B] a new code was developed which also use similar techniques 
to handle the axis problem. This article is particularly relevant in our con- 
text because they use essentially the same equations and the same gauge 
conditions as presented in section [2j However, the mass formula (l39i) for this 
gauge was not known at that time and hence it was not used in that article. 

Further studied were performed in [29J [31J, where new evolutions schemes 
are proposed and the general case with twist is considered. In [32] the reg- 
ularization procedure at the axis was performed for more general evolutions 
schemes. In [30] important issues concerning the maximal-isothermal gauge 
are studied, in particular a new choice of variable is proposed in order to 
guarantee the uniqueness of solutions of the relevant equations. I will come 
back to this in the next section. 

In a very recent series of articles, the critical collapse of gravitational 
waves is studied [33] [M]. This code use a different gauge choice, namely 
harmonic coordinates. Finally, we mention [20] where the linear equations in 
the maximal-isothermal gauge were studied and boundary conditions com- 
patible with the mass formula were presented. 
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3.2 Analytical results 



The first analytical application of the mass formula is the proof of the fol- 
lowing geometrical inequality 



where J is the total angular momentum of the spacetime (see [H] [T7] for 
details). Extension of this theorem, which remove technical assumptions and 
also simplify considerable its proof were presented in [8] [9] . The remarkable 
aspect of this proof if that extreme Kerr black hole appears as the minimizer 
of the mass integral. This suggests stability properties of these black holes. 
In particular, it suggests that it is perhaps easier to prove stability near 
an extreme black holes than a non-extreme one in axial symmetry (see the 
essay [IS]). There exist important extensions to inequality (jUl) which include 
charges [UJ [TO] . 

The second geometrical inequality proved with a suitable quasi-local for- 
mulation of the mass formula is the following 



where A is the area and J is the angular momentum of the black hole horizon. 
See [18] , pQ [21] for details. In the stationary case with matter and charge, it 
has been proved in [21] . Numerical evidences for the validity of this inequality 
has been presented in [25J. 

Inequality (J4"2"]l ensures that the Christodoulou quasi local-mass of the 
black hole is a monotonically increasing quantity in the evolution (see the 
discussion in [TS]). 




Both geometrical inequalities are proven on the initial data, that is, the 
proofs do not make use of the evolution equations. In this gauge the equations 
reduce to a coupled hyperbolic-elliptic system which is formally singular at 
the axis. Due to the rather peculiar properties of the system, the local in time 
existence has proved to resist analysis by standard methods. To analyze the 
principal part of the equations, which may represent the main source of the 
difficulties, we study linear perturbation around the flat Minkowski solution 
in this gauge. In [22] we solved this linearized system explicitly in terms of 
integral transformations in a remarkable simple form. This representation is 
well suited to obtain useful estimates to apply in the non-linear case. 

4 Open Problems 

There exists two relevant open model problems for axially symmetric space- 
times: 




(41) 



8tt\J\ < A 



(42) 
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1. Black hole formation and critical collapse of axially symmetric gravi- 
tational waves (Numerical). 

2. The stability of the Kerr black hole in axial symmetry (Mathematical). 

These problems are expected to be difficult, they constitute the long term 
motivation of this study. In the following I will briefly discuss them in order 
to propose simpler open problems that can contribute to the resolution of 1 
and 2 and at the same time are feasible to solve in the near future. 

Most of the numerical work mentioned in section l3TTl is concerned with the 
vacuum equations and with the critical collapse of the gravitational waves. 
As we pointed out in the introduction, there exists other interesting physical 
systems in axial symmetry. For example, the collapse of an axially symmetric 
rotating star studied in [26] . However, if we restrict ourselves to the vacuum 
equations (as we do in this article) then probably the most relevant open 
problem to study in numerical axial symmetry is the critical collapse. And 
the reason is that this problem can not be handled, with present computers, 
without symmetry assumptions due to the high resolution needed. The sub- 
ject has already a long history (see the review article [5]), and despite many 
efforts the problem still lacks definitive answer (see the recent article |34|). 

Although there exist many possible gauge choices in axial symmetry, 
the mass formula suggests the maximal-isothermal gauge as privileged one. 
Also, constrained evolutions schemes appear to be better than free evolu- 
tions schemes (see the discussion in [30]). In the following we focus on the 
maximal-isothermal gauge in its constrained formulation. In the article [30J 
it was clearly pointed out one of the main problems of these equations. There 
are many elliptic equations and it is not clear a priori that they will always 
admit a unique solution along the evolution. A new set of variables was pro- 
posed in that article to ensure that the equations have the correct behavior 
at the linearized level to ensure uniqueness of solution. However, it is still 
not clear if the non-linear equations (notably, the Hamiltonian constraint) 
is always solvable. It is an open problem to see if these equations (or some 
variant of them) are well posed. I believe that this is an important problem 
that stress the interplay between analytical and numerical studies. Perhaps 
the resolution of the well-posedness question will lead us to select (or even 
discover) the correct evolution scheme. 

Another question which is important in numerical evolution is the issue 
of outer boundary conditions on a finite domain. In the article [20J a set of 
boundary condition at the linear level were proposed. This boundary condi- 
tions have the property that the energy appears to leave the domain, at least 
for the class of initial data studied in this article. It is important to note that 
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the mass formula provides a well-defined local measure of the gravitational 
energy. These boundary conditions have not yet been implemented for the 
complete equations. 

The following two problems summarize the discussion above. 

l.a. Well-posedness for the maximal isothermal gauge. IstheHamil- 
tonian constraint always solvable? Does the gauge breaks down in the 
strong field regime? 

l.b. Implementation of the radiation boundary conditions for the 
full equations. How much they improve the evolution scheme? It is 
possible to prove (may be with some extra assumptions) that for such 
boundary conditions the energy always leaves the domain? 

As we have seen, problem 1 is more numerically oriented, since it ap- 
pears to be very difficult to understand critical phenomena with the current 
analytical tools. In contrast problem 2 is more analytical. In recent years 
there have been important new developments in the problem of the linear 
stability of the Kerr black hole from the analytical point of view. Most of 
this work is concerned with scalar waves on a fixed Kerr black hole back- 
ground. We refer to the review articles [12] and [13] and reference therein. 
These works are not restricted to axial symmetry. The expectation is, of 
course, that axially symmetric spacetimes are simpler. In particular, axial 
symmetry provides extra geometrical structure that can give new insight into 
the equations. As we mentioned in section f3T2l much of the current evidences 
that this is the case come from the geometrical inequalities (jUJ) and fj42|) . 
These inequalities have been proved under some restricted conditions. In the 
case of the quasi-local inequality ( H2"j) . the only restriction is that the initial 
data are assumed to be maximal. It is important to remove this restriction 
or to find counter examples. Numerics can play a relevant role in the search 
of possible counter examples. Also, it would be very interesting to generalize 
this inequality to include the charge. This is important not only because the 
charge is a relevant parameter of a black hole but also because the charge is 
well defined without any symmetry assumption (this is not the case of the 
quasi- local angular momentum). The natural question is whether a similar 
inequality holds without axial symmetry. This would involve an appropriate 
definition of quasi-local angular momentum. Another interesting problem is 
to study what happens in the case of equality in (H2"|) . A rigidity property 
has been proved in [2T] that characterize the local geometry of the surface 
in that case. However, it is very likely that this kind of surface can only ap- 
pear as asymptotic limit as the extreme Kerr initial data. A related question 
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is to study small deformation of the extreme Kerr initial data, as the ones 
constructed in [19J. Do these initial conditions admit apparent horizons if 
the perturbation is small enough? This question is relevant because if the 
answer is negative then these initial conditions have the chance to provide 
complete Cauchy surfaces that are outside the black hole region, as in the 
extreme Kerr black hole. 

Regarding the global inequality OH]), the most relevant open problems is 
to prove it for multiple black holes and also to remove the maximal condition. 
The multiple black holes case is related with the uniqueness of the Kerr black 
hole among stationary black holes with disconnected horizon components. 

All the previous problems are concerned with the initial conditions. The 
maximal-isothermal gauge will be useful to answer questions regarding the 
black hole stability if the mass formula can be used in some way to control the 
evolution. A natural first problem would be to recover the non-linear stability 
of Minkowski [7] in this gauge. The expectation is that the mass formula will 
provide a simpler (and different) kind of approach to this problem; although, 
of course, always restricted to axial symmetry. A second problem it is the 
study of axially symmetric perturbations of a black hole in this gauge. The 
expectation is that the mass formula can help to prove linear stability under 
axially symmetric perturbations of the Kerr black hole. The advantage of 
this problem in comparison with the non-linear stability of Minkowski in this 
gauge is that it is probably easier to deal with and also it will prove something 
new, since the linear stability of the Kerr black hole is still an open question. 

We summarize this discussion in the following list of open problems in 
relation to problem 2. 

2. a. Geometrical inequalities. Remove the maximal assumption. For 
the quasi-local inequality (jl2"l) include the charge and study the possible 
generalization of this inequality without any symmetry assumptions. 
Study the case of the equality: it is possible to have such a surface 
in an asymptotically flat initial data? Study the existence of minimal 
surfaces or apparent horizon for small perturbation of extreme Kerr 
initial data. For the global inequality (j4ip . prove it for multiples black 
holes. 

2.b. Linear stability of the Kerr black hole in axial symmetry. Use 

the mass conservation to control the linear evolution. 
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